Stability of quasi-two-dimensional Bose-Einstein condensates 
with dominant dipole-dipole interactions 
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We consider quasi- two-dimensional atomic/molecular Bose-Einstein condensates with both con- 
tact and dipole-dipole interactions. It is shown that, as a consequence of the dimensional reduction, 
and within mean-field theory, the condensates do not develop unstable excitation spectra, even when 
the dipole-dipole interaction completely dominates the contact interaction. 
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Dilute Bose-Einstein condensates with long-range in- 
teractions offer promising opportunities to explore the 
potentially strong correlations induced by the interac- 
tion, which go beyond the comparatively weak correla- 
tions induced by the local contact interaction pseudopo- 
tential conventionally sufficient to describe most atomic 
condensates. The recent realization of a concrete physi- 
cal system in which atomic magnetic dipoles play a sig- 
nificant role Q, is a first experimental step towards the 
exploration of dipolar condensate physics 0, S IS @ • 
To investigate dipole- dominated physics, greater poten- 
tial than by magnetic dipoles is offered by heteronuclear 
molecules with an electric dipole moment 0, 0| , due to 
the fact that the electric dipole interaction strength is 
larger than the magnetic one by a factor of about 10 4 
for typical atomic/molecular dipole moments of 1 Bohr 
magneton or 1 Debye, respectively. 

An interesting property of three-dimensional (3D) con- 
densates is the existence of a "roton" minimum, which 
stems from the fact that the Fourier transform of the 
dipole-dipole interaction, occurring in the Bogoliubov 
mean-field spectrum, assumes negati ve values in certain 
directions in momentum space jilllCt fll| . This property 
however also causes a problem of 3D dipolar condensates, 
as they become dynamically unstable for dipole-dipole 
interactions dominating the contact interaction, at den- 
sities corresponding to the Thomas-Fermi limit of large 
condensates. The roton minimum quickly deepens with 
increasing density and/or dipole coupling to hit the zero 
of energy, and a further increase is impossible because 
the system develops imaginary excitation energies. 

In the present study an analytical proof is given that, 
in contrast to 3D condensates, quasi-two-dimensional 
(quasi-2D) dipolar condensates are stable, even when the 
dipole coupling completely dominates the contact cou- 
pling. The experimental realization of purely dipolar 
condensates is thus possible only in the quasi-2D regime. 
The physical reason for the stabilization is that dimen- 
sional reduction entails that the mutual dipole-dipole in- 
teraction of the atoms/molecules, with dipoles oriented 
perpendicular to the plane in the strongly confined di- 
rection, can effectively sample an exponentially smaller 



region in configuration space where the interaction as- 
sumes negative values. It is, furthermore, shown that 
a strongly correlated regime of purely dipolar conden- 
sates, where large quantum depletion of the condensate 
occurs and the mean-field description in terms of a single 
condensate wave function breaks down, takes place in a 
crossover regime from quasi-2D to 3D. For an electric- 
dipole-dominated dilute gas, the densities corresponding 
to this crossover regime to 3D turn out to be quite small, 
even within the strongly enlarged stability window of 
strongly anisotropic pancake-shaped condensates. 

The mean-field description of the system we start with 
is based on a nonlocal 3D Gross-Pitaevskii equation for 
the order parameter 'J (we put H = m = 1, where m is 
the mass of the atoms/molecules) 0,0] : 



^- + ^y+g m \^(r,t)\ 2 (1) 



d 3 r'V dd (r-r')\*(r',t)\ 2 \H>(r,t). 



The trapping in the plane is assumed to be negligible 
as compared with the strong confinement, of harmonic 
trapping frequency lo z , in the z direction. The term in the 
second line contains the dipole-dipole interaction of the 
atoms/molecules, with all dipoles oriented by an applied 
field along the strongly confining z direction: 



V dd (r) 



3g d 1 - 3z 2 /|rf 
47r Irl 3 



(2) 



The dipole coupling reads g d — /io^m/3 for magnetic and 
g d = dg/3eo for electric dipoles. In the units which we 
employ, the g d are having dimensions of length like the 3D 
contact interaction g^D = Aira s , for which the length scale 
is set by the s-wave scattering length a s . We note that 
the value of the contact interaction coupling strength 
generally depends on the value of g d (in an effective- 
dimension dependent manner), because the long-range 
dipole-dipole interaction affects the short-range scatter- 
ing processes . For the statements made in the present 
paper, the actual ratio 33D / g d and the absolute value of 
the (unscreened) g d will be relevant. 
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For a quasi-2D Bose-Einstein condensate, the motion 
of the atoms/molecules is by definition restricted to zero 
point oscillations in a harmonic oscillator potential. We 
then take as a general ansatz for the density 



pW = |*(r)| 5 



f 



: exp 



z 



n(x,y), (3) 



where the density in the plane, n(x,y), is normalized 
to the total number of particles N. We treat d z as a 
parameter minimizing the Gross-Pitaevskii ground state 
energy 0, El- Assuming homogeneous density in the 
plane, the equation determining d z reads u> 2 = d~ + 
(#3D + 2g d )nNd~ 3 /V^n. If the right-hand side is dom- 
inated by the first (kinetic energy) term, d z equals the 
harmonic oscillator length the system is quasi- 

2D, and the above Gaussian gives the density profile ex- 
actly. Defining a parameter a = uj z d 2 , we have a = 1 
if the system is quasi-2D and a 3> 1 deep into the 3D 
regime, where a becomes interaction dependent. 

We now calculate the total dipole-dipole energy given 
that the density profile in z-direction is prescribed by the 
above Gaussian. In accordance with Eq. JTJ, the Gross- 
Pitaevskii energy functional of the dipole-dipole interac- 
tion generally reads 



2 
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2 (2tt) 3 



d 3 r / d 3 r'p{r)V dd (r - r')p(r') 



d 6 kp(k)V dd (k)p{-k), 



(4) 



where the second line employs a convolution to Fourier 
space. The Fourier transform of the dipole-dipole interac- 
tion © takes the form V dd (k) = g d [3k 2 J ' {k 2 x +k 2 y +k 2 z )~V\- 
using J2Jl, we obtain the Fourier transform of the density 
p{k) = exp [—jk 2 dl] n(k x ,k y ). Integrating over the k z 
direction, the effective dipole-dipole energy is then given 
by a two-dimensional integral in (k x , k y ) space 
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where the complementary error function Erfc [z] = 1 — 
erf (z) = 1 - (2/Vtt) Q exp(-t 2 )dt and k = (k 2 x + k 2 ) 1 / 2 ; 
h{k x ,k y ) is the Fourier transform of the 2D density. 
Employing the same procedure of integrating out the 
Gaussian in the z direction, for the contact interaction 
part in the Gross-Pitaevskii energy functional, yields the 
well-known result for the 2D effective coupling ^2D = 
2\/27ra s /d z (valid in the limit that the 3D s-wave scat- 
tering length a s <C d z fl^l. 

From the relation 10 for the dipole-dipole contribu- 
tion in the interaction energy, it follows that the Fourier 
transform of the total, contact plus dipole, interaction 




FIG. 1: Fourier transform of the total interaction potential in 
units of A I p(0)y / 7vd z , as a function of £ = kd z . The values of 
R = \^ft[2 (solid line), R = § v / V/2 (dotted line) R = 
(dipole interaction dominating; dashed line). For R = |^/7r/2 
(i.e., when gd = g3o), the Fourier transform approaches zero 
at infinite momentum. 



potential assumes the form 
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where we made use of the w-function, related to Erfc(z) 
by w(z) = exp [z 2 ]Erfc(z) [H Q. The dimcnsionlcss 
wavenumber £ = kd Zl and the two dimcnsionlcss param- 
eters occurring in ttot > (0; using the central 3D density 
p(0), are defined to be 



A = 



p(0)y/^d 2 z g d 
R 



R 



1 + 33D/2.9d' 



(7) 



The value of the parameter R ranges from R = if 
9d/93B -> 0, to R = iJtt]2 for g d /g 3 B -> oo. The 
second term in the curly brackets in Eq. J^J rapidly de- 
creases as a function of £ and approaches a constant, 
which is due to the fact that w((,/\/2) — > 1/2/77 £ _1 
for ( -> 00, cf. the plot of V?J?(() in Fig.Q] For 
small C, the quasi-2D Fourier transform behaves like 
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C + ^§C 2 + 0(C 3 ) 
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serve that V^(Q possesses a well-defined value at the 
origin ( = k = 0, as opposed to the 3D Fourier transform 
of the dipole-dipole interaction potential. 

From the Fourier transform of the interaction |JBJ, 
we conclude that the squared Bogoliubov spectrum 
[Isl flfij |. for excitations confined to the plane, oj 2 = 
p(0)^d z V 2 ^(k)k 2 + k 4 /A is, in units of 1/d*, given by 



V2. 



4 ' 



(8) 



The main observation of the present analysis is that, as 
opposed to the 3D case (or, in an exacerbated manner, 
the quasi- ID case 0), the Bogoliubov spectrum (jHJ does 
not necessarily become unstable if the dipole interaction 
coupling exceeds the contact interaction coupling. As- 
suming 930/ g d — * 0, i.e. R — > ■u , tt/2, the above squared 
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FIG. 2: Squared Bogoliubov spectrum in units of l/d*, 
Eq. ©, where R = (gzT>/g d 0) for all curves. Count- 

ing from bottom to top at small kd z the values of A in Eq. Q 
are: A = A c /10, A = A min , A = ,4 C , and A = 1.2A C . For 
yl > j4 m i n = 1.249, the spectrum displays a roton minimum, 
and becomes unstable for A > A c = 3.446. 

spectrum can assume negative values, and hence the ex- 
citation energies become imaginary, when A exceeds the 
critical value A c — 3.446, cf. Fig.[3 where we show 
a sequence of four spectra for different A at constant 
R = spK]2. 

The critical value of A c — 3.446 corresponds to a crit- 
ical dipole coupling given by (gd) c — 2.436/ p{Q)d 2 z . Us- 
ing numbers appropriate for atomic magnetic moments 
Mm = N m (i B , we find that gd — 5.4 x 10~ 6 /imMA^ 
(M is the mass of the atoms or molecules in units 
of the atomic mass unit = 1.66 x 10~ 27 kg). In the 
case of electric dipoles, with moment d e = N e Debye, 
we find g d = 6.2 x 10~ 2 pmMN 2 . The critical cou- 
pling (gd) c for dominant dipole-dipole interactions then 
translates into a critical central 3D density p c (0) = 
4.5 x lO^cm^iV^a-^pTr x kHz] and p c (0) = 3.8 x 
10 12 cm _3 iV e _2 Q; _1 Ci; 2 [27r x kHz] in the case of magnetic 
and electric dipoles, respectively. 

The quantity ^/np(0)gd/R, a measure of the energy 
per particle (the chemical potential), is equal to A/d 2 . 
For the system to be quasi-2D, we therefore need A <C 
to z d 2 = a = 1. Furthermore, the Bogoliubov spectrum 
in Fig.El does not possess any points where de 2 /d( 2 =0 
for values A < A min = 1.249 (when R = y/n/2). That 
is, a "roton" minimum cannot develop within the regime 
of quasi-2D. Thus we conclude that a quasi-2D purely 
dipolar system of bosons is always stable, and that the 
instability of the condensate takes place in the crossover 
region to 3D. By contrast, deep inside the 3D regime, 
.93D "C gd implies collapse of the condensate for any rea- 
sonable value of the density, and in the Thomas-Fermi 
limit the condensate will be unstable for any gd slightly 
exceeding g w 0. 

The critical value of A C (R), at which the gas becomes 
unstable, obtained by numerically finding the large mo- 
mentum zeros of the Bogoliubov spectrum JHJ, exponen- 
tially increases for smaller R (increasing gsv/gd), and 
diverges at R — |^/ 7r /2 (i-e., for g^c = gd), where 



the Fourier transform V 2 ^ becomes positive everywhere. 
Comparing the critical values of A thus obtained to those 
from a full 3D solution of the Bogoliubov-de Gcnncs 
equations [9J, the present approach reproduces the criti- 
cal A sufficiently accurate in the dipole-dominated case, 
for which the instability takes place in the crossover 
regime from quasi-2D to 3D. Deep into the 3D regime, 
the critical value A c calculated from JHJ) is strongly over- 
estimated, mainly because it is exponentially sensitive on 
the exact form of the spectrum. 

At nonzero temperature T, 2D Bose-Einstein conden- 
sates do not exist in the homogeneous case |l7|]. while 
trapping in the plane enables the existence of a conden- 
sate also at finite T O n the other hand, in the 
presently discussed T = case, it is a well-established 
fact that even without trapping, Bose-Einstein conden- 
sation occurs in two spatial dimensions, see, e.g., (lij. We 
next turn to a discussion of the zero temperature value of 
the number density of excitations above the condensate, 
the so-called quantum depletion. To this end, we use a 
mode expansion for the annihilation operators \k of the 
original bosons in terms of the Bogoliubov quasiparticle 
operators afc,a^ |20| : 




The above form of the Bogoliubov transformation results, 
after inversion, in the usual phonon quasiparticle oper- 
ators at low momenta, and gives = at fc — > cx), 
i.e., the quasiparticles and the bare bosons become, as 
required, identical at large momenta. 

The quantum depletion density at zero temperature is 
calculated by evaluating the expectation of y^X m the 
quasiparticle vacuum defined by afc|vac) = 0: 

The instability generally happens because in-plane ex- 
citations become of imaginary frequency 0; above the 
calculated critical value of A = A c , the in-plane excita- 
tions have vanishing energy, and the quantum depletion 
diverges. Using the in-plane momentum integral above 
to calculate the depletion, with the effective dispersion 
relation JSJl, assumes that the corrections due to the ne- 
glected out-of-plane Bogoliubov excitations, character- 
ized by transverse quantum numbers, are small. This is 
justified because the dominant contribution to quantum 
depletion is, for a dilute Bose gas, from large momentum 
excitations with (approximately) vanishing energy. Out- 
of-plane excitations have large energies at the relevant 
in-plane momenta of order l/d z ; they thus do not con- 
tribute significantly to the depletion. In Fig.|3 we show 
the result for the quantum depletion in the case of dipole- 
interaction dominated condensates, for which we can ex- 
pect the in-plane spectrum (JSJ) to be sufficiently accurate 
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FIG. 3: Quantum depletion density from Eq. 110L as a func- 
tion of A = p(0)y^ndlg d /R, for R = y/ n/2 (g d /gw — > oo). 



Hall states of purely dipolar Fermi gases were studied in 
[2^ . The interplay of quasi-two-dimensionality and rota- 
tion may generally yield interesting new physics [3] . It 
will, furthermore, be of interest to determine the change 
in the system's stability region when the dipoles are no 
longer locked in one direction, and thus to investigate the 
influence of spin waves on the stability of dipolar quan- 
tum gases of bosons. 

I thank C. Zimmermann, N. Schopohl, L. Santos, C. 
Iniotakis, and N. R. Cooper for helpful discussions. 



up to A ~ A c . At the critical value A = A c ( y / ir/2) ~ 3.4 
(dotted vertical line in Fig. 01, the condensate depletion 
diverges, and the mean-field condensate will yield to a 
new quantum phase. 

In conclusion, we have shown that dipolar quasi-2D 
Bose-Einstein condensates are extremely stable systems 
as compared to their 3D counterparts. This offers the 
potential of approaching, starting from the mean-field 
physics of condensates, a strongly correlated regime of 
dilute atomic/molecular gases with long-range interac- 
tions. A conceivable experimental procedure is to start 
from a quasi-2D dipolar condensate, and, by decreasing 
u) z , to enter the crossover regime to 3D, where the quan- 
tum depletion of the condensate becomes large. If the 
dipole-dipole interaction dominates, i.e., in the limit R — > 
V^72, the quantity A/V2 = p(0)d 2 z g d = 2Ap(0)/p c (0), 
where /o c (0) is the critical 3D density for Bogoliubov ex- 
citations above the mean-field condensate to be stable, 
discussed in the above. Hence p(0)d 2 z gd measures the di- 
luteness of the system, that is the proximity to the critical 
region. Numerical values are 



p(0)d 2 z9d 



5.5xio-^ y( y"rJ 

oj z [2tt x kHz] 

Ar2 p(0)[10 12 cm- 3 ] 

0.63aW e =- , 

lo z [2ir x kHz] 



(11) 



in the magnetic and electric cases, respectively. Accord- 
ing to the first equation, magnetically dipolar gases will 
allow for the realization of a purely dipolar condensate, 
such that the depletion still remains small. On the other 
hand, it is evident from the second line of Eq. 111(1 . that 
for electric dipoles the densities at which strong deple- 
tion of the condensate sets in are rather small, even for 
strongly increased axial trapping. 

Among further lines of research offered by the stability 
potential of quasi-2D dipolar condensates are the phe- 
nomena expected when they are set in rotation. When 
both contact and dipole interaction are present, various 
of these phenomena have been explored in |2l| ; quantum 
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